Astronomy & Astrophysics manuscript no. aa ' SpMicFlare2006 ' 20final ' figsrenamed © ESO 2008 

February 2, 2008 
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ABSTRACT 



We present an analysis of contemporaneous photospheric, chromospheric and coronal structures on the highly active K-dwarf star BO 
Mic (Speedy Mic). We concentrate on a moderate flare that we localize in the stellar atmosphere and study its energetics, size and 
thermal behavior. The analysis is based on strictly simultaneous X-ray, UV- and optical observations carried out by XMM-Newton 
and the VLT. We use Doppler imaging and related methods for the localization of features. Based on X-ray spectroscopy we study 
the the coronal plasma in and outside the flare. The flare emits in X-rays and UV, but is not detected in white light; it is located at 
intermediate latitude between an extended spot group and the weakly spotted pole. We estimate its height to be below 0.4 stellar radii, 
making it clearly distinct in longitude and height from the prominences found about two stellar radii above the surface. In contrast to 
BO Mic's photospheric brightness, neither its chromospheric nor its X-ray emission show a pronounced rotational modulation. 

Key words, stars: activity - chromospheres - coronae - flare - late-type - stars: imaging - individual: BO Mic 
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1. Introduction 

The level of stellar activity is related to the efficiency of a 
star's dynamo, operated by an interplay of convection and ro- 
tation, with f aster rotation resulting in h igher activity for a given 
type of star dSchriiver & Zwaanl I2000). Extreme activity phe- 
nomena are observed on the ultrafast rotating, highly active, 
young K-dwarf star BO Mic (K2V, P ro , = 0.380 + 0.004 days, 
vsin/s; 135km/s, d - 44.5 pc), nicknamed "Speedy Mic". 
BO Mic produced the largest X-ray flare observed during the 
Rosat all-sky survey, with a peak flu x of » 9 • 10 31 erg/s in the 
Rosat PSPC passband (0.1 - 2.4 keV. iKiirsterll 1995b . Follow-up 
observations of BO Mic showed its very high projected rota- 
tional velocity and short rotation period. Subsequent Doppler 
imaging studies revealed a richly spotted photosphere and re- 
configur ations of active re gions during o nly a few ste llar ro- 
tations (I Wolter et al J 120051 = WSW 2005. lBarnesll2005h . mak- 
ing strictly simultaneous and coherent observations mandatory. 
We observed BO Mic during two consecutive nights on 2006 
October 13/14 and 14/15 simultaneously with XMM-Newton and 
Vlt/Uves at ESO Paranal. The observations were scheduled to 
maximize simultaneous and continuous visibility of the target 
by both instruments. 

In this letter we report the first results of our campaign, pre- 
senting a coherent "snapshot" of BO Mic's active regions that 
concentrates on the localization and analysis of a moderate flare. 

2. Observations and data analysis 

Our Vlt/Uves observations lasted two half nights, JDo + 0.49 
to 0.70 and JD + 1.48 to 1.69 with JD = 2454022.0. The 
XMM-Newton observations exceeded these time spans, lasting 



* Based on observations obtained at the ESO VLT Obs. No. 078.D- 
0865(A) mid XMM-Newton Obs. Id. 0400460301 and 0400460401 



from from JD + 0.48 to 0.79 and from JD + 1.47 to 1.73. 
In total we obtained about 50 ks of X-ray data and 32 expo- 
sures of optical photometry from the optical monitor (OM) on- 
board XMM-Newton. Our 142 Uves spectra completely cover 
two semi-rotations of BO Mic with some phase overlap. 

Uves was operated in dichroic mode, covering the spectral 
range 3260 to 9460 A with a spectral resolution of A/AA » 
40 000. Uves exposures lasted between 200 and 250 s, separated 
by the CCD readout time of 10 s. The spectra were reduced using 
the package Reduce described in lPiskunov & Valentil d2002l) . We 
carried out an approximate flux calibration using synthetic spec- 
tra from Eiuichndritii] dH). 

All XMM-Newton instruments obtained useful data of 
BO Mic; here we only present data gathered by the EPIC 
(European Photon Imaging Camera) and the OM (optical mon- 
itor). The OM was operated with the visual grism in the wave- 
length range 3000-7000 A. Due to the brightness of BO Mic, 
the OM grism spectra are affected by coincidence loss, and we 
could not make use of either the full spectral resolution or of 
the absolute flux calibration, thus we had to restrict the analy- 
sis to relative spectrophotometry. All XMM-Newton data were 
reduced using the S cience Analysis System (SAS) version 7.0 
(lLoiseau et al.l l2.006h with standard selection and filtering cri- 
teria. X-ray light curves are background subtracted, based on 
nearby source-free regions. Spectral analysi s of the EPIC data 
was carried out with XSPEC version 11.3 (lArnaudlll996l) ; for 
spectral fitting purposes we used a multi-temperature model, as- 
suming the emission spectrum of a collisionally i onized, opti- 
cally thin gas as calculated with the APEC code dSmith et al.l 
2001) with abundances modeled once for the complete data set. 
To investigate spectral changes we then derived temperatures 
and volume emission measures (EM= J n e nndV) of the individ- 
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ual plasma components and calculated X-ray luminosities from 
the best fit models. 

3. BO Mie s activity spatially resolved 

3.1. Doppler images 

The photospheric spots were reconstructed from the Vlt spec- 
troscopy using our Doppler imaging (DI) algorithm Cldi 
(WSW 2005) on a grid resolving 2020 surface elements. We 
adopted an inclination angle of 70° for the stellar rotation 
axis, v sin = 134 km/s a s determined by WSW 2005 in agree- 
ment with Barnes (2005), and a rotation period of 0.380 days 
dCutispoto et al.in~997l) . For DI, the spectra were added in pairs, 
resulting in 72 spectra with a typical S/N ratio of 400 at 6000 A. 
The spectra cover 1 . 1 stellar rotations with a homogeneous phase 
sampling and a gap of two rotations between the observations 
of the two hemispheres. The line profiles used for the DI were 
obtained by a spectrum deconvolution of the wavelength range 
6390-6440 A (see WSW 2005 for details). 
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Fig. 1. Atmospheric structures of BO Mic reconstructed by 
Doppler imaging and related methods. All maps show the same 
features for the given rotation phases. Photospheric spot cover- 
age is rendered as black, dark and light orange areas, represent- 
ing 100%, 67% and 33% spot coverage, respectively. Equatorial 
arrows mark approximate longitudes of prominences observed in 
Ca K and H a . The blue circle indicates the chromospheric emis- 
sion region of the flare. Grid lines mark latitudes and longitudes 
in 30° steps. 

In Fig. [T] we show a reconstruction of BO Mic's active re- 
gi ons. The spots on the surface show a largely non-polar (similar 
to iBarnesI d2005l) . although at lower latitude), azimuthally quite 
asymmetric distribution, with spots both in the northern and 
southern hemispheres. The large-scale reliability of our Doppler 
image, at least in longitudinal spot distribution, is illustrated in 
Fig. [2] where we compare the observed light curve with that 
computed from our DI in the 6000-6500 A range. The light 
curve supports the reality of the southern spots found in the DI. 
While their appearance in the DI depends on the adopted spot 
brightness, the reconstructions including southern spots yield a 
significantly better fit to the OM-light curve. 
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Fig. 2. BO Mic's photospheric emission as measured by 
XMM-Newton 's optical monitor (black) and computed from our 
Doppler image (blue). The left and right sections of the black 
curve render the first (lower phase scale) and second (upper 
scale) observed semi-rotation, respectively. The emission is in- 
tegrated from 6000 to 6500 A , i.e. centered between Johnson- V 
and R. The estimated error is +0.003 mag. 
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Fig. 3. BO Mic's coronal and chromospheric emission. Purple 
crosses show flux and lcr-errors of the soft X-ray emission (0.2- 
5.0 keV) as registered by XMM-Newton 's EPIC-pn detector. The 
black and red curves give the Ca K and H,, equivalent width, 
respectively. Arrows mark the phases of sub-observer passage 
for the prominences shown in Figs.[T]and [5] The horizontal lines 
"A-G" indicate bin intervals for the data points of Fig. [4] Phase 
annotations are the same as in Fig. [2] 



Concerning the non-periodic components of the OM-light 
curve, we note that the peaks at phases 1 .25 and 4.25 differ in 
brightness. This reflects a change in the spot pattern during the 
three rotations between them, as supported by changes in the 
corresponding photospheric line profiles. 

3.2. The chromospheric and coronal perspective 
In Fig. |3] we plot the temporal evolution of the chromospheric 
and coronal activity in terms of Ca K and H ff equivalent width 
and soft X-ray emission in the 0.2-5 keV energy band (crosses). 
Clearly, neither varies smoothly between the two nights. The X- 
ray light curve shows flares during both nights; the second night 
seems to consist of one coherent event with several bumps in the 
decaying X-ray light curve, suggesting reheating events. 

To verify this interpretation we investigated the evolution of 
the X-ray emitting plasma in terms of emission measure and 
temperature. We divided the X-ray data into seven segments ("A- 
G" in Fig. [3]) and performed a spectral analysis (cf. Sect. [2]). The 
results are shown in Fig. [4] the temperature decreases with de- 
creasing emission measure during the intervals A-D, which is 
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consistent with the decay phase of a flare. During interval E the 
temperature rises, indicating additional energy input. The tem- 
peratures of 20-30 MK are typical for moderate flares. 
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Fig. 4. Evolution of emission measure and temperature for the 
hot plasma component during the flaring, with lcr errors as de- 
termind by the X-ray spectral fit. EM 5 is used as a proxy for the 
plasma density. The solid line shows a linear fit to points B-D, 
corresponding to the decay of the first flare event. Further energy 
release is evidenced by the deviation of points E-G from the de- 
cay path. The time bins of points "A-G" are given in Fig. [3] 

Outside the flaring phases the X-ray light curve is relatively 
flat and apparently dominated by a stable quasi-quiescent level; 
it varies only by 20% peak to peak. As visible in Fig. [3] this 
modulation has a local maximum centered on phase (p « 1 .55, so 
it may be correlated with the photospheric brightness minimum 
at (f> « 1.65 (Fig.|2]l. However, based on one observed rotation, 
and in the presence of flaring, this remains speculative. 
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Fig. 5. Evolution of the Ca K emission: Ca K line core pro- 
files, processed by unsharp masking are rendered as a function 
of wavelength (units of v sin i) and rotation phase. Absorption 
dips in the profile appear blue while emission bumps are shown 
red. Plusses indicate sine-fits according to Eq. [T] They are used 
to localize absorbing prominences, their longitudes are given to 
the right. The fit to the flare-induced chromospheric emission is 
marked by diamonds. See text for discussion. 



3.3. Localizing prominences and the flare 

In Fig.|5]we provide a detailed view of the evolution of the Ca K 
emission line profile during our observations. Specifically, we 
show the time series of line profiles processed by an unsharp- 
masking filter that enhances fast changes in the core profile. The 
unsharp masking is performed by dividing the profile time series 
by a temporally smoothed copy of itself. For Fig. [5] wavelengths 
were transformed to velocity units, absorption dips in the profile 
are shown in blue while emission features are shown red. Several 
absorption features clearly move through the Ca K line profiles. 
In order to localize the corresponding prominences in longitude 
we fit pronounced deformations migrating redward through the 
core profiles with sine functions depending on rotation phase <p 
through 

v rad(0) — v sin i ■ r/R„ cosinjT. — 8) sin(ip — 2n<p) (1) 

which describes the radial velocity of an atmospheric feature lo- 
cated at the radius r (divided by the stellar radius R„), the lat- 
itude 8 and longitude tp. Essentially, the phase of v rac i = de- 
termines the longitude while the slope ^ determines the height 
of a prominence when 8 is adopted as the subobserver latitude. 
In this way we obtain the approximate longitudes and heights of 
six prominences marked in Fig. [5] These prominence parame- 
ters agree with those obtained by the same procedure from the 
H a core profiles. Their radii range between 2.5 and 3.5 R* with 
an uncertainty of +0.5 R t , i.e. they are about 2 R, above the sur- 
face. In Figs. Q] and [3] the approximate longitudes and phases 



of the structures causing absorption transients in the Hydrogen 
Balmer and CaH and K lines are marked by arrows. Such 
fast-moving, often periodic line profile deformations can be at- 
tributed to "prominences", i.e. clouds of relatively cool material 
magnetically kept in co-rotation h igh above the stellar surface , 
and already observed on BO Mic dCollier Cameron et al.| [2003: 
iDunstone et alj|2006l) . Given their rapid movement through the 
line profile, these prominences cannot be located close to the 
stellar surface. Only one pronounced deformation of the Ca K 
line core moves through the profile compatible with an emission 
feature close to the stellar surface. Shown in Fig. [5] it passes the 
profile center at phase (f> = 4.0 + 0.02. Using Eq. Q] to trace the 
emission peak through the line profile, we locate its emitting re- 
gion on the surface at a latitude of 8 — 56 + 10° and a longitude 
of if — 359 + 10°, marked by a circle in Fig.Q] 

It should be noted that this "flare plage" could also be located 
at the same latitude south of the equator, since for each rotation 
phase the radial velocity is the same there as for the northern 
location. The same effect causes a north-s outh-ambiguity of fea - 
tures found in all Doppler images (e.g. Pisku nov et alj fl990). 
However, in that case, the weak emission rise in H ff and Ca K at 
phase <p ss 4.3 (see Fig. [3} could not originate from the plage lo- 
cation because south of the equator this latitude is no longer vis- 
ible. Thus, this feature is very likely located in the north. Since 
the enhanced Ca K and H ff emissions clearly coincide with the 
X-ray flare we are confident that the location marked in Fig.[T]is 
indeed the chromospheric footpoint of the X-ray flare. 
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4. The flare: evolution and energetics 

Our contemporaneous optical and X-ray data allow us to further 
characterize the flaring plasma. Unfortunately, our observations 
did not cover the impulsive rise of the chromospheric and coro- 
nal flare emission and we see no flare signatures in the photo- 
metric lightcurves. However, a short increase prior to the peak 
of the X-ray emission is visible in Fig. [3] suggesting that our ob- 
servations cover the very end of the rise phase of the flare and all 
of its decay. The decay to the "quiescent" level takes about 0.45 
and 0.29 rotation phases or 4.1 and 2.6 hours for the soft X-ray 
and Ca K emission, respectively. The H a emission decays even 
faster, however, it is strongly influenced by the prominences. 

Assuming that the flare oc curs in a single p lasma loop and 
using the method described in iReale et all d2004l) . we estimate 
the flare loop size from the EM-T diagram shown in Fig. [4] 
Fitting the decay of the first heating event, the derived slope 
yields an upper limit of I ;S 0.4 R„ for the loop half length. This 
upper limit even remains valid if the flare occurred in an arcade 
of several loops, because in this case each individual loop would 
have to be smaller. While considerably larger than solar coro- 
nal loops, this is much closer to the star than the prominences 
located at h « 2R* (cf. Sect. 13 . 3 b - Reheating occurs between 
rotation phases 4.2 and 4.3, indicating its connection to the sud- 
den slight increase in X-rays seen in Fig. [3] at phase 4.25. This 
is possibly related to the increase in H a and Ca K emission ob- 
served about 0.05 rotations (« 30 min.) later. The total emitted 
energy of the flare observed in soft X-rays, ^ 1.3- 10 34 erg in the 
0.2-10 keV band, appears to originate predominantly, but not 
exclusively, from a single energy release. The peak luminosity 
of the flare in soft X-rays of w 1.4 • 10 30 erg/s (6 x 10 29 erg/s in 
the 1-8 A GOES-band) is considerably larger than the maximum 
power in Ca K and H„ of ~ 2.0 • 10 28 erg/s and « 1 .0 • 10 29 erg/s, 
respectively. This energetic mismatch of coronal and chromo- 
spheric emission would be unusual for a solar flare where these 
emissions - with considerable variations - tend to be comp arable 
(lFoukalll2004 Uohns-Krull euMi99% lEmslie et alJl2005t) . 

5. Summary and discussion 

We observed the highly active K-star BO Mic during one ro- 
tation with completely simultaneous coverage from soft X-ray 
to optical wavelengths. The observations were carried out in 
two blocks separated by two stellar rotations with all observed 
indicators of photospheric, chromospheric and coronal activity 
showing pronounced variations. 

The Doppler image (DI) shows only a few spots near the 
visible pole. Most spots are asymmetrically distributed at mid- 
latitudes, leading to a roughly sinusoidally modulated photo- 
spheric light curve with a peak-to-peak amplitude of w 0.1 mag. 
There are additional small non-periodic components in the op- 
tical light curve; in agreement with results of previous observa- 
tions of BO Mic (WSW 2005), this indicates some reconfigura- 
tions of spots on rather short timescales of a few rotations. 

The chromospheric light curves (Hq, and Ca II) and coro- 
nal light curve (soft X-ray emission) are only weakly correlated 
with the photospheric modulation. The predominant variations 
of chromospheric emission are absorption transients in the H a 
and Ca K line cores, which are caused by prominences about 
two stellar radii above the surface; their distribu tion is qualita- 
tively similar to those analyzed for BO Mic by Dunstone et al. 
(2006). The longitudes of the prominences show no obvious 
connection to the photospheric spots. Whether the "quiescent" 
chromospheric emission is more intensive on the photospheri- 
cally darker hemisphere, as suggested by our 2002 observations 



dWolter & Schmittll2005h . cannot be reliably determined due to 
the strong prominence absorptions found here. Also, there is no 
clear evidence for X-ray absorption by these prominences, set- 
ting a strict upper limit to their column density. 

The most conspicuous event during our observations was a 
flare lasting about 4 hours; the flare is not registered in our broad 
band photometry, but shows up clearly in the H,,, Ca II and X-ray 
light curves. By tracing the radial velocity of its chromospheric 
emission in the Ca K line we localize the flare site in the stellar 
atmosphere. Somewhat surprisingly, the flare does not appear to 
be connected to the main concentration of activity in terms of 
dark spots; rather, it is located at the fringes of this region. As 
discussed in Sect. 13.31 this "flare plage" (which could also be a 
post-flare loop) could also be situated at the same latitude south 
of the equator, where the DI is only poorly defined due to poor 
visibility close to the stellar limb. At any rate, the "flare plage" is 
not associated with any distinguished feature of the photospheric 
spot pattern. Either such a spot feature is largely hidden by the 
southern limb or, more likely, it does not exist. 

We determine the overall energetics and thermal evolution 
of the flare by simultaneously studying the soft X-ray emission. 
With a peak soft X-ray power of « 10 30 erg/s, this flare is a hun- 
dred times more energetic than a la rge solar flare in terms o f 
peak flux and total radiated energy dGolub & PasachoffHl997l) . 
however, it is still at least two orders of magnitude weaker than 
the largest flares observed on BO Mic. Assuming « 10 30 erg/s 
as an upper limit to the power radiated in the optical continuum, 
this would correspond to ~ 10~ 3 of BO Mies total luminosity 
(WSW 2005), i.e. below the sensitivity of our photometry and 
consistent with the photometric non-detection. 

In addition to the large flare, BO Mic's soft X-ray emission 
exhibits another, smaller flaring event and possibly a slow mod- 
ulation with * 20% peak to peak amplitude. The slow modu- 
lation may be correlated with the observed photospheric spot 
distribution; however, due to the flares and only one observed 
rotation, this remains inconclusive. Also, no optical spectra are 
available for the first flare, therefore it could not be localized. 
However, both X-ray flares appeared at similar rotation phases 
(</> * 2.0 and 4.0), i.e., the first flare was situated on the same 
hemisphere as the second or sufficiently close to one of the poles. 
Qualitatively, our X- ray lightcurve is similar to that observed by 
iHussain et all d2007l) for the ultrafast rotating K-star AB Dor. 
Also, our value for the coronal flare loop height of $ 0.4/?* is 
compatible with their height estimate of coronal structures, once 
more confirming the similarity of AB Dor and BO Mic. 

To our knowledge, the flare on BO Mic is the first X-ray 
flare that could be localized on a Doppler image. Interestingly, 
the flare occurs on a rather inconspicuous portion of the atmo- 
sphere as judged from the image. Extended observations of flar- 
ing events on a given star could reveal the spatial distribution 
of flare emerging sites. Physically, one expects flaring events to 
be associated with magnetic null lines. Since our UVES spectra 
are non-polarimetric, we cannot assess the location of magnetic 
null lines on BO Mic's surface. Clearly, a Zeeman Doppler im- 
age would be helpful to prove the association of stellar flare sites 
with magnetic null lines as expected from the solar analogy. 
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